On the existence of a bound tetraneutron 
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Following recent work in which events which may correspond to a bound tetraneutron ( 4 n) were 
observed, it is pointed out that from the theoretical perspective the two-body nucleon-nucleon force 
cannot by itself bind four neutrons, even if it can bind a dineutron. A very strong phenomenological 
four-nucleon (4N) force is needed in order to bind the tetraneutron. Such a 4N force, if it existed, 
would bind 4 He by about 100 MeV. Alternative experiments such as ( 8 He, 4 n) are proposed to search 
for the tetraneutron. 



PACS numbers: 21.45,+v, 21.10.Jx, 25.60.Je, 27.10.+h, 
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In a recently reported experiment Jj], ||] events were ob- 
served that exhibit the characteristics of a multineutron 
cluster liberated in the breakup of 14 Be, most probably 
in the channel 10 Be+ 4 n. The lifetime of order 100 ns 
or longer suggested by this measurement, would indicate 
that the tetraneutron is particle stable. If confirmed, 
the existence of a bound tetraneutron is of great im- 
portance as it would challenge our understanding of nu- 
clear few-body systems and nucleon-nucleon (NN) inter- 
actions. Here I would like to make several observations 
in this context. 

I 

As discussed in Refs |], ^ the breakup 14 Be — > 
10 Be+ 4 n represents an attractive channel to search for 
a tetraneutron. The nucleus 10 Be is strongly bound and 
the four neutron separation energy for 14 Be is only about 
5 MeV with respect to the 14 Be — ► 10 Be + 4n breakup. 
It is known that in an external potential well neutrons 
can form drops PL 0. In the attractive field of 10 Be the 
last four neutrons could form a tetraneutron which might 
be liberated in the breakup. 

In earlier experiments, the tetraneutron was 
searched using heavy-ion transfer reactions such as 

and double ex- 
These reactions 
should be strongly suppressed by spectroscopic reasons 
and a negative outcome from them could be easily 
anticipated. 
II 

No proper ab-initio four-body calculations of the tetra- 
neutron with realistic two-body and three-body NN 
forces are known to the author. However, several the- 
oretical calculations of the tetraneutron are available. 

(i) In Ref. || the tetraneutron, studied in the trans- 
lationally invariant, symmetrized oscillator basis using 
the Afiuj model space, was found to be unbound by 18.5 
MeV. However, the oscillator basis is not appropriate for 
the description of unbound systems. 

(ii) No bound tetraneutron was found in Ref. JT(i| within 
the angular potential functions method with semirealistic 
NN interactions. No search for the four-body resonance 
state was carried out there. 

{iii) No bound tetraneutron was found in Ref. p"l| ] within 
the stochastic variational method on a correlated Gaus- 
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sian basis for a range of simple effective NN potentials. 
No search for the four-body resonance state was carried 
out there. 

(iv) A search of the four-body resonance in the lowest or- 
der of the hyperspherical functions method (HSFM) gave 
a null result |l|, ||]. 

(v) The energy behaviour of the eigen phases, studied in 
Ref. @ within the HSFM using the K max = 6 model 
space, led the authors to the conclusion that the tetra- 
neutron may exist as a resonance in the four-body con- 
tinuum at an energy of about 1 — 3 MeV. However, such 
a conclusion is not convincing because no convergency 
of the eigen phases with an increase of the hyperangular 
momentum has been achieved in these calculations. Be- 
sides, a clear indication of the resonance has been seen 
with only one of the NN potentials used in the calcu- 
lations, namely with the Volkov effective NN force VI 
p^f . Volkov effective NN interactions reproduce the ex- 
perimental binding energy of another four-body system, 
4 He. However, their singlet even and triplet even com- 
ponents are equal, and therefore with these potentials a 
singlet dineutron has exactly the same binding energy as 
the deuteron. In the particular case of VI, the dineu- 
tron is bound by 0.547 MeV. Therefore, VI cannot be 
used in calculations of the multineutron systems. An- 
other NN potential, used in Ref. |L4j , namely that of Re- 
ichstein and Tang (RT) |ll|, reproduces the n-p triplet 
and p-p singlet scattering lengths and does not bind a 
dineutron. With this potential the energy derivatives of 
the eigen phases monotonically decrease with increasing 
energy. This means that resonances in the four neutron 
system are absent, at least within the model space con- 
sidered. 

To understand whether the RT potential can produce 
any resonance or bound state if the model space of the 
HSFM is increased, I have calculated the hyperradial po- 
tentials for the tetraneutron up to K max = 16 using 
the technique developed in Ref. fTij] . The RT poten- 
tial has only central even components Vio and Voi and 
the odd potentials are usually obtained from them us- 
ing an arbitrary parameter u so that Voo = (u — l)Voi 
and Vii = (u — l)Vio- The present calculations have 
been performed with u = 1 which corresponds to no- 
interaction in the odd partial waves. The diagonalised 
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hyperradial potentials Vdiag(p), calculated for different 
model spaces K max , are plotted in Fig.l as a function of 
hyperradius p. Such potentials can be used to obtain a 
quick approximate solution for the binding energy in the 
extreme adiabatic approximation of the HSFM [081 . One 
can see that the Vdiag(p) have almost converged. They 
are purely repulsive, monotonically decreasing with the 
hyperradius, and do not show any sign of local attractive 
pockets. Therefore, the RT potential can neither bind 
the tctraneutron nor produce any resonances. To get a 
bound tetraneutron, the value of u should be increased 
up to u — 2.3. However, with this value, the isotopes 4 H 
and 5 He must become bound with respect to the t + n 
and 4 He + n thresholds by about 6.5 MeV and 32 MeV 
respectively. In reality, these nuclei are unbound by 3 
and 0.9 MeV. 
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FIG. 1: Diagonalised hyperradial potentials Vdi ag (p) of the 
tetraneutron calculated with Kmax ~ 2, 4, 16 with two 
different central NN potentials. 

For comparison, the calculations of the hyperradial 
potentials have been performed with Volkov NN poten- 
tial VI as well. The diagonalised hyperradial potentials 
Vdiag(p), shown in Fig. 2, reveal local attractive pock- 
ets for K max > 6. These pockets become negative for 
K m ax > 12, but they are too shallow to form a bound 
state with respect to the four-body decay. Although with 
K m ax = 16 convergence has not yet been reached, the 
general trend seen in Fig. 2, suggests that it is unlikely, 
that a tetraneutron, bound with respect to the 2 n + 2 n 
decay, may exist. To get a bound tetraneutron, the Ma- 
jorana parameter m should be changed from its standard 
value of 0.6 to m = —0.2. Such a change provides E( 4 n) 
= —1.2 MeV and does not influence the binding energy 
of 4 He, however, it binds 4 H and 5 He with respect to 
the t + n and 4 He + n thresholds by 7 MeV 25 MeV 
respectively. 

And finally, calculations with the realistic two-body 
isospin-conserving central NN interaction Argonne i>4 [p~S| ] 
have been performed as well and the corresponding hy- 
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FIG. 2: The same as in Fig.l for the Volkov potential VI. 



perradial potentials Vdiag(p) are shown in Fig.l. These 
purely repulsive potentials have almost converged and 
their behaviour excludes any possibility of either a bound 
or a resonance state. 
Ill 

The results of the theoretical calculations suggest that 
the tetraneutron must be unbound and most likely should 
not exist as a resonance, at least when only two-body 
central forces are considered. It is remarkable that even 
when an effective NN potential binds the dineutron (the 
case of VI), it still cannot bind two dineutrons, although 
a resonance in such a system should exist. In order for 
the two-body central force to be able to bind the tetra- 
neutron, a huge unphysical attraction should be intro- 
duced in the triplet odd potential. Such an attraction 
would strongly overbind 4 H, 5 He and other known A > 4 
nuclei. The reason for necessity of huge attraction lies 
in the relative number of nucleon pairs in even and odd 
states allowed by the Pauli principle. For the tetraneu- 
tron with L = and S = the probabilities P to find a 
pair of nucleon in the singlet even (01) and triplet odd 
(11) states are equal, fni = Pu = |- For the tetra- 
neutron with L = 1 and S = 1, these probabilities are 
■Poi = 3 an d Pu — §■ The NN force is not sufficiently 
attractive in the singlet even state, and is repulsive in 
the triplet odd state, which leads to the absence of the 
tetraneutron. For comparison, for the L — 0, S = state 
of the closed shell nucleus 4 He Pxo = Poi = \ and for 
the L = 1, S = 1 state of 4 H P w = P 01 = Pu = \. The 
binding energies of 4 He and 4 H are —28.3 MeV and —5.2 
MeV respectively. It is clear that the binding of a nucleus 
is strongly correlated with the probabilities to find a pair 
of nucleons in even states, especially in the triplet even 
state (10) where the n — p bound state exists. 

IV 

The fact that, despite the theoretical predictions made 
with central two-body forces, six events 14 Be — > 10 Be + 
4 n have possibly been observed, means that either non- 
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FIG. 3: The DWBA a-transfer cross sections of the 
d( 8 He, 4 n) 6 Li reaction as a function of the tetraneutron bind- 
ing energy. 



central interactions and/or three- nucleon (3N) forces, 
and/or a strong four- nucleon (4N) force may bind the 
tetraneutron. If not, then a different interpretation of 
the six events from Refs pi |2| should be looked for. 

The first possibility should be rejected, because the 
calculations of Ref. jl(| indeed included the spin-orbit 
and tensor forces. The second possibility does not look 
convincing either because the contribution of the 3N po- 
tential to the binding energies of the neutron drops n 
and 8 n, calculated with Green's function Monte Carlo 
Methods in Ref. [^0| , is about 1 to 5 MeV which is small 
compared to 4 He and to other A = 7 and A — 8 nuclei. 
A similar contribution from the 3N force would be too 
small to bind the tetraneutron. Therefore, if what has 
been seen in the 14 Be experiment was indeed the tetra- 
neutron, it could indicate at presence of a 4N force. 

At present, the results of the ab-initio calculations of 
the 3 < A < 8 nuclei suggest that a 4N force is not 
needed to fit the observed binding energies of light nuclei 
at the 1% level @. Therefore, either the 4N contribu- 
tions are smaller than 1% for these nuclei or parts of the 
3N forces are mocking up their effects |2(J . On the other 
hand, the latest models of the 3N force used to calcu- 
lated the binding energies of these nuclei have not yet 
been tested in the description of polarization observables 
of the low-energy Nd scattering where the contribution 
of the 3N force is very important [S. Besides, the 
latest numerically accurate calculations employing a 3N 
force do not reproduce the proton analyzing power for 
the p- 3 He scattering p5|. The simultaneous fit of the 
Nd and p — 3 He polarization observables and of binding 
energies of the lightest nuclei may leave room for a 4N 
force. 

In order to get an idea of what strength is needed for 
the 4N force, the HSFM calculations of this paper have 
been repeated with the RT potential using a 4N force 
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FIG. 4: The same as Fig. 3, but for the 12 C( 8 ii'e, 4 n) 16 
reaction. 



simulated by the potential V^m{p) = Woe~ ap . The val- 
ues a = 0.7, 1.2 and 1.5 /to -1 , used in calculations, were 
the same as the range of the phenomenological term of 
the 3N spin-orbit force introduced in Ref. [^lj . The tetra- 
neutron becomes bound if the corresponding values Wq 
are equal to —410, —1460 and —2530 MeV respectively. 
These values are two orders of magnitude larger as com- 
pared to the values —1, —10, and —20 MeV obtained 
for the 3N spin-orbit force [|lj. For the Volkov poten- 
tial VI, these strengths are similar: Wq = —320, —1565 
and -2900 MeV. If the same 4N force existed in 4 Hc, 
the binding energy of 4 He, calculated with VI, would be 
-88, -82 and -134 MeV respectively. Therefore, the 
4N force, if it exists, should be strongly T-dependent. In 
principle, such a strong T — 2 4N force could be noti- 
cable in A — Z > 4 nuclei, especially near the neutron 
edge of stability, where the number of the T = 2 four- 
nucleon states increases. However, HSFM calculations of 
the neutron-rich helium isotopes 6,8 He, performed with 
VI [[[7), do not indicate a need for a strong 4N force. 

To get a real understanding of the role of the 4N force, 
ab-initio calculations of the tetraneutron with realistic 
2N and 3N forces must be performed. The HSFM is the 
best method for such calculations because it is well suited 
to systems without bound subsystems, which is the case 
for the tetraneutron, and because it can also treat a four- 
body continuum. 

V 

Since the discovery of the tetraneutron would most cer- 
tainly require a revision of modern realistic models of the 
three-nucleon interaction, it is extremely important to 
confirm or reject the results reported in Ref. There- 
fore, independent searches for the tetraneutron in alter- 
native reactions should continue. Among such reactions, 
a transfer reactions A( 8 He, 4 n)A + a deserve special at- 
tention. Using a light target and detecting a product 
nucleus A + a, one can study not only a hypothetical 
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tetraneutron bound state but also the four-neutron con- 
tinuum as well. 

The spectroscopic factor S for the ( 8 -ffe| 4 -ffe <8> 4 n) 
overlap can be estimated using the translation invariant 
shell model p4]. If the probability amplitudes of the 
\L = 0, S = 0) and \L — 1,S — 1) configurations are 
cto and ai in 8 He and /3o and /3i in n respectively, then 

S = (y^aoA) + ^/fai/Si) ■ 

The choice of the target A must be determined by 
kinematic and spectroscopic considerations. It should be 
light enough for residual nucleus A + a to come out of 
the target and its wave function should overlap strongly 
with the A + a wave function. The obvious targets of 
this type are the deuteron, triton, 3 He and 4 He. The list 
of other light targets with large spectroscopic factors S a 
is given in Ref . pf| . 

As an example, in Figs. 3 and 4 the DWBA a-transfer 
cross sections of the d( s He, 4 n) 6 Li and 12 G{ 8 He, 4 n) 16 
reactions are shown as a function of the four-body 
binding energy of the hypothetical tetraneutron. The 
energies of these reactions have been chosen based on 
availability of optical potentials in the entrance and 
exit channels, which have been taken from |2(| and 
f27f . In addition, the energy choice for the 12 C target 
was also motivated by experimental measurements of 
the similar reaction 12 C( 8 Li, a) 16 N which revealed very 
large cross sections at small angles, about 25 mb/sr 



The latest theoretical value S a = 0.87 Q has 
been used for 6 Li and S a for 8 He was taken to be |. 
According to Figs. 3 and 4, the d( 8 He, 4 n) 6 Li reaction 
is better suited for the tetraneutron search because its 
cross sections are not strongly influenced by unknown 
tetraneutron binding energy. The cross sections and 8 He 
beam intensities make this experiment feasible in the 
very near future. In contrast, the cross sections of the 
12 C( 8 -ffe, 4 n) 16 reaction decrease dramatically with 
increasing the binding of the tetraneutron and become 
small if the tetraneutron is bound. 

Summarizing, due to the large probability for a pair 
of neutrons to be in the triplet odd state, the two-body 
nucleon-nucleon force cannot by itself bind four neutrons, 
even if it can bind a dineutron. A very strong phe- 
nomenological T-dependent four-nucleon force is needed 
in order to bind the tetraneutron. Alternative experi- 
ments such as ( 8 He, 4 n) are proposed to search for the 
tetraneutron. 
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